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ABSTRACT 


The  linear  stability  of  an  incompressible  laminar  flow  in  the  blending  boundary  layer 
between  the  boundary  layer  in  a  90'’  streamwise  corner  and  a  Blasius  boundary  layer  well 
away  from  the  corner  is  examined  using  a  locally  parallel  flow  approximation.  It  is  shown 
that  the  influence  of  the  outer  boundary  conditions  associated  with  oblique  modes  of  dis¬ 
turbances  which  are  anti-symmetric  about  the  bi.sector  plane  have  a  profound  effect  on  the 
stability  of  the  flow.  As  a  result,  in  good  agreement  with  observation,  the  critical  streamwise 
Reynolds  number,  associated  with  a  spanwise  location  is  significantly  reduced  as  the  corner 
is  approached,  being  Rct  =  60  approxitnately  for  spanwise  distance  of  c*  =  from  the 

corner  compared  with  Rcr  =  322  approximately  for  c*  =  20j;*/?“’,  where  measures  down¬ 
stream  distance  from  the  leading  edges.  At  R  =  600,  the  growth  rate  of  the  most  amplified 
mode  of  disturbance  at  the  former  location  is  over  six  times  greater  than  that  at  the  latter; 
the  corresponding  wave  angle  at  the  two  locations  is  respectively  44°  and  5°,  approximately. 
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1.  Introduction 

The  secondary  flow  which  arises  in  the  streaniwise  corners  of  a  wind  tunnel  has  important 
bearing  on  the  size  of  the  available  test  section,  in  experiments  involving  laminar  flow  in  a 
two-dimensional  boundary  layer,  for  example,  the  role  of  the  side  walls  of  the  tunnel  in  the 
observed  transition  is  not  always  clear.  Further  the  spread  of  the  secondary  flow  from  the 
junction  of  a  wing  with  the  fuselage  of  an  aircraft  has  implications  for  laminar  flow  control 
on  wings. 

A  knowledge  of  the  nature  of  the  flow  in  a  streaniwise  corner  and  its  stability  character¬ 
istics  is  therefore  of  considerable  practical  importance.  A  model  flow  which  has  been  studied 
is  that  in  which  the  corner  is  formed  by  two  semi-infinite  plane  surfaces  inclined  at  an  angle, 
usually  a  right  angle,  to  each  other  with  the  free  stream  flow  along  the  corner  (figure  1). 
The  basic  flow  has  a  similarity  solution  which  is  symmetric  about  the  bisector  plane  (see  Ru- 
bin(1966),  for  example)  such  that  well  away  from  the  corner  it  matches  the  two-dimensional 
Blasius  boundary  layer  flow  along  the  plane  surfaces  while  in  the  corner  it  is  governed  by 
fully  three-di?iiensional  equations  of  motion.  In  between  these  regions,  Rubin(1966)  identi¬ 
fied  a  ‘blending  boundary  layer’  where  the  influence  of  the  corner  decays  algebraically  away 
from  the  corner.  Further.  Rubin  showed  that  to  leading  order  the  spanwise  velocity  profile 
in  this  region  exhibits  a  reverse  flow  in  the  boundary  layer. 

Zamir(I98I)  estimated  from  experimental  evidence  that  in  the  absence  of  pressure  gra¬ 
dient,  transition  in  a  right  angled  streaniwise  corner  occurs  at  a  value  of  the  square  root  of 
the  Reynolds  number,  based  on  the  streaniwise  distance,  of  R  =  100  (henceforth  R  will  be 
referred  to  as  the  Reynolds  number).  Uidike  the  experimental  results  of  El-Gamal  and  Bar- 
clay(1978),  Zamir  and  Young  (1979)  found  breakdown  of  similarity  in  the  velocity  profiles. 
Zamir(1981)  attributed  this  as  a  development  of  flow  instability. 

Lakin  and  Hussaini  (1984)  derived  the  equations  governing  small  perturbations  of  the 
basic  flow  in  the  ‘blending’  boundary  layer  identified  by  Rubin  using  a  locally  parallel  flow 
approximation.  They  considered  heuristic  solutions  of  these  equations  based  on  the  critical 
layer  in  the  boundary  layer. 

In  this  paper,  we  re-consider  the  stability  of  the  blending  boundary  layer  for  a  right 
angled  streaniwise  corner,  using  a  locally  parallel  flow  approximation  as  done  by  Lakin  and 
Hussaini  and  show  that  the  instability  is  dominantly  driven  by  the  outer  boundary  conditions 
cissociated  with  disturbances  which  are,  unlike  the  basic  flow,  anti-symmetric  about  the 
bisector  plane.  The  spanwise  cross-flow  turns  out  to  be  of  too  insignificant  a  magnitude  to 
account  for  the  observed  early  transition  in  the  flow,  unlike  in  the  case  of  flow  jiast  a  swept 
wing,  for  example. 


In  section  2,  the  equations  of  motion  and  the  boundary  conditions  are  given.  The  equa¬ 
tions  obtained  are  exactly  those  given  by  Lakin  and  Hussaini  (1984)  although  it  has  not 
been  necessary  to  make  any  a  priori  assumptions  about  the  nature  of  the  disturbances.  The 
equation  governing  the  normal  component  of  the  perturbed  velocity  satisfy  a  modified  Orr- 
Sommerfeld  equation.  Since  the  basic  flow  is  symmetric  about  the  bisector  plane,  the  outer 
boundary  conditions  need  to  be  applied  on  the  bisector  plane  rather  than  at  ‘infinity’.  This 
gives  rise  to  the  possibility  of  modes  of  disturbances  which  are  symmetric  or  anti-symmetric 
about  the  bisector  plane.  It  is  found  that  the  anti-symmetric  modes  are  the  most  unstable. 

The  method  of  orthonormalization  is  used  to  determine  the  eigenvalues  of  the  stability 
problem  and  is  briefly  described  in  section  3. 

The  results  of  the  investigation  are  presented  in  section  4.  As  expected,  it  is  found 
that  the  critical  Reynolds  number,  decreases  with  decrease  in  spanwise  distance  from 
the  corner,  it  being  given  by  Rcr  =  54  for  a  spanwise  distance  of  compared 

with  R  =  322  for  2*  =  20/?"^ 3:*;  the  growth  rate  of  the  the  most  unstable  disturbance  at 
R  =  600  is  over  six  times  greater  at  the  former  location  than  that  at  the  the  latter  while 
the  corresponding  wave-angle  is  over  nine  times  greater.  Further,  it  is  found  that  in  the 
proximity  of  the  corner,  the  three-dimensional  modes  of  disturbance  are  unstable  at  a  lower 
Reynolds  number  than  that  for  two-dimensional  ones. 

2.  Basic  Equations  and  boundary  conditions 

We  define  Cartesian  co-ordinates  {x*,y*,z*)  as  shown  in  figure  1,  with  the  two  semi¬ 
infinite  rigid  planes  intersecting  at  right  angle  in  the  line  Ox".  A  viscous  incompressible 
fluid  of  uniform  density  flows  along  the  corner  with  the  undisturbed  velocity  away  from  the 
planes  given  by  (t/*,0, 0)  where  U‘  is  a  constant.  The  undisturbed  flow  is  symmetric  about 
the  bisector  plane  y*  =  2*  so  that  it  is  only  necessary  to  consider  the  flow  in  the  region 
between  this  bisector  plane  and  the  rigid  plane  y*  =  0,  say,  with  appropriate  boundary 
conditions  on  the  two  planes.  The  fluid  velocity  satisfies  the  exact  non-dimensional  equations 
(u*  =  U*u{x,y,  z),  t*  =  x*  =  lx,p*  =  p*U*^p  where  I  denotes  the  downstream 

distance  of  the  location  under  consideration): 

— -fu.Vu  =  -Vp-I-— (1) 

and 

V.u  =  0,  (2) 
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where  R  =  ([/*l f is  the  Reynolds  number.  The  boundary  conditions  to  be  satisfied  are 


u(a:,0,2:,t)  =  0,  u(x,oo,oo,t)  =  (1,0,0), 

together  with  appropriate  boundary  condition  on  the  bisector  plane  x  = 
We  consider  a  small  perturbation  about  a  basic  .steady  flow,  so  that 

u  =  u(x)  4-  u(x,  <),  p  —  p{x)  +  p(x,  t). 


(3) 


(4) 


On  substituting  (4)  into  (1)  and  respectively  equating  zeroth  order  and  first  order  terms  in 
u  and  p  to  zero  we  have 


u.Vu  +  Vp  =  — 

— +;u.Vu  +  u.Vu  =  —Vp-i-—V^u 
V.u  =  0,  V.u  =  0 

with,  writing  u  =  {U,  V,  W), 

u(x,0,2)  =  0,  u(x,oo,oo)  =  (1,0,0) 

au  ac  „  „ 

_  =  — ,  v  =  Wor,!,  =  . 

and 


(5) 

(6) 
(7) 


(8) 

(9) 


u(x,0,2,t)  =  0,  u(x,oo,oo)  =  0 

together  with  appropriate  boundary  conditions  on  u&iy  —  z. 

In  region  S  <  z  «  x,  where  S  is  the  non-dimensional  thickness  of  a  Blasius  boundary 
layer,  S  =  4.9x’/^/?“’,  Rubin  (1966)  has  shown  that  for  a  steady  flow,  symmetric  about  the 
plane  y  =  z,  (5)  has  an  approximate  similarity  solution. 


!=(/(, )  +  0(fi-'-), R-'(u!(.))  +  0(-)))  +  0[R-‘) 


(10) 


X  X'* 

where  a  prime  denotes  differentiation  with  respect  to  the  ocpendent  variable  r/  =  yR/y/x 
and  V  =  lipf'  —  f),w  =  aog{Ti)  and  /  and  g  satisfy 


/  +  2//  =  0;  /(O)  =  /  (0)  =  0;  /'((X))  =  1 
9”  +  \if9)'  =  fif(O)  =  0;(7(oo)  =  -1 


(11) 


(12) 


with  oo  =  .86039.  Here  the  streamwise  pressure  gradient  is  o{R  ').  Thus  the  approximate 
basic  flow  corresponds  to  a  Blasius  boundary  layer  flow  with  an  imposed  cross-flow,  R~^w. 
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associated  with  the  normal  component  of  velocity  in  the  boundary  layer  on  the  other  rigid 
plane. 

We  consider  the  location  x  =  l(i’*  =  /),  ~  =  R~^zo  (~o  >  b),  and  introduce  the  concept 
of  locally  parallel  flow.  Without  loss  of  generality,  we  consider  a  disturbance  centered  on 
(l,y,f2“'2o)  on  an  appropriate  scale  of  0( /?“').  We  write 


X  =  l  +  y  —  R  ^  7]\  0{  R  z  —  R  ^  {zq z\),  t  —  R  ^  t\  (Id) 

and  note  that  //  =  r/i  +  V  =  ftV]  etc.  Thus  locally,  in  the  vicinity  of 

{l,y,R-^zo), 

u  =  {f'{m)^R-'viin)^R-'Myu))  +  0{R-^).  (14) 


If  we  write  u  —  (n,  l>,u;),  then  on  substituting  (14)  and  (15)  into  (6)  and  (7),  we  obtain  after 
some  manipulation, 


r,  ^  1  i)  R  I  II  d  ,  _ 

du  dv  dw 

u - 1”  15 - 1"  15 —  “ 

(Jx\  u1)\  (JZ\ 

- ,  -//  dv  I  I  dv  . 

where  C  =  ~  ^  's  the  normal  component  of  perturbed  vorticity  and 


Equations  (15)-(17)  have  a  solution  of  the  form 


(15) 


(16) 

(17) 


(18) 


(19) 


where  v\  and  Ci  satisfy  the  de-coupled  equations; 

[[D^  —  —  [vD +  v  ^- iR{af  ■]- —  cic)){D^  —  -\- iR{oif  + )]t>]  =  0  (20) 

/t  H 

[*«(«/'  +  ^w-ac)  +  vD  -(D^-  A2)]C,  +  iRidf"  -  ^w')]v,  =  0.  (21) 

K  rC 

Here  D  =  d/dzx  and  4-  /i^.  Equation  (20)  corresponds  to  the  Orr-Sommerfeld 

equation  for  a  basic  flow  f/  =  iZ  +  it  is  usual  to  drop  the  terms  involving  nandn  ,  but 

these  are  retained  here  consistent  with  retaining  R~^w  in  the  basic  flow.  Retaining  the.se 
terms  makes  minor  changes  in  the  stability  characteristics.  Equation  (20)  is  the  same  as 
given  by  Lakin  and  Hussaini  (1984)  although  in  deriving  it,  it  has  not  been  neces.sary  to 
make  any  a  priori  assumption  about  the  nature  of  the  disturbance. 
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Boundary  conditions 


The  approximate  basic  flow  in  (20)  is  independent  of  the  spanwise  location  c  and  it  is 
tempting  to  obtain  (as  is  usually  done  in  the  case  of  the  Orr-Sommerfeld  problem  associated 
with  the  Blasius  boundary  layer)  the  outer  boundary  conditions  from  (9)  and  from  the 
limiting  form  as  ;/i  — >■  oo  of  (20)and  apply  these  at  a  large  value  of  r/j ;  the  limiting  form  of  the 
fourth  order  equation  (20)  has  two  exponentially  growing  solutions  and  two  exponentially 
decaying  solutions,  one  solution  in  each  category  being  viscous  and  the  other  inviscid  in 
character.  Care,  however,  is  required  in  the  present  case  since  the  basic  flow  is  regarded 
as  being  symmetric  about  the  bisector  plane  y  —  z  so  that  the  outer  boundary  condition 
must  strictly  be  applied  at  this  plane,  instead  of  at  some  large  value  of  r/^  as  the  restriction 
4.9  <  Zo  <<  R  imply  that  the  normal  to  the  plane  y  =  0  at  the  spanwise  location  r  = 
meets  the  bisector  plane  at  rji  =  rj^out  =  -o  i  where  4.9  <  i]\out  «  R-  It  is  through  this 
requirement  that  the  dependence  of  the  stability  characteristics  on  the  spanwise  location 
enters  into  the  problem. 

Although  the  basic  flow  is  symmetric  about  the  plane  y  =  z,  no  such  restriction  applies 
to  the  perturbations.  Accordingly,  the  perturbation  velocity  u  =  (u,t>,tv)  may  satisfy  one 
of  four  possible  basic  sets  of  boundary  conditions  on  y  =  z:  (i)§^  =  §z,  ^  =  u.b(ii)  = 
V  =  —w,  (iii)  u  =  0,  V  =  w,  (iv)  u  =  0,  v  =  —w.  Conditions  (i)  and  (ii)  imply  that  u  is 
symmetric  while  v  ~  tu  is  respectively  anti-symmetric  and  symmetric  m  y  =  z.  Conditions 
(iii)  and  (iv)  imply  that  u  is  anti-symmetric  while  v  —  w  \s  respectively  anti-symmetric  and 
symmetric  in  y  =  z.  It  is  found  that  the  most  unstable  mode  of  disturbance  corresponds 
to  condition  (iii).  We  therefore  consider  only  this  case  in  detail  here;  the  other  cases  may 
be  considered  similarly  without  any  difficulty.  In  view  of  the  continuity  equation,  it  follows 
from  condition  (iii)  that  we  must  have 

{D  -b  il3)v  -  0  on  t/,  =  cq  (-2) 

Next,  we  note  that  in  spite  of  the  restriction  on  its  range,  7]\out  still  lies  outside  the  boundary 
layer  so  that  we  may  consider  the  limiting  form  of  (20)  as  t/i  — +  >>  1  in  determining  the 
second  outer  boundary  condition.  Hence, 

{D'^  -  A^)(D  -  x)v  =  0  on  T/i  =  ro  >>  1  (23) 

where  X  =  |(u(oo)  —  (itJ^(cx))  +  4(A^ -}- z7?(o'  + /yu;(oo)//?)))’^^).  The  choice  (23)  nu’ansthatof 
the  four  solutions  to  the  limiting  form  of  (20)  described  above,  we  only  discard  the  viscous 
exponentially  growing  solution,  the  growth  rate  of  this  solution  being  much  greater  than  the 


corresponding  inviscid  solution.  The  conditions  at  the  surface  t/i  =  0  are 


n,(0)  =  Dn,(0)  =  0 


(24) 


3.  Numerical  method 


The  linear  equation  (20)  together  with  the  boundary  conditions  (22)-(24)  is  solved  using 
the  technique  of  orthonormalization  (see  for  example,  Davey  (1973))  which  allows  the  use 
of  standard  shooting  methods.  We  first  express  (20)  as  a  system  of  first  order  ordinary 
differential  equations; 

it  =  \,..A)  (25) 

where 


<f>2  =  <j>\,  <f>3  = 

4>a  =  (f>3^  (26) 

and  the  non-zero  elements  of  the  matrix  a,j  are: 

<Il2  =  ^23  =  034  =  1?  O21  = 

041  =  —iR{af  "  +  f R) 

043  =  +  iR{af  +^w/R  —  ac)  +  v,  044  =  v.  (27) 

The  boundary  conditions  (22)-(24)  take  the  form: 


4>2{^o)  +  —  0 

^4(20)  -  X<t>3i-o)  =  0 

<^i(0)  =  <^2(0)  =  0  (28) 

Equations  (25)  are  to  be  integrated  from  rji  =  zq  for  various  values  of  zq  to  r/i  =  0.  The 
steady  flow,  characterized  by  /',  v  and  w,  is  evaluated  at  2M  points  by  integrating  (14)-(17) 
using  fourth  order  Runge-Kutta  integration. 

Corresponding  to  the  two  outer  boundary  conditions  given  in  (28),  two  independent 
solutions  satisfying  these  conditions  and  the  four  equations  (25)  can  be  chosen.  We  denote 
these  by  vectors  {j  =  1,2)  where 
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with  starting  values  at  7/1  =  zq  given  by: 


^<’>  =  [1, -7/^,0, of 

^<^)  =  [0,0,l,xf  (29) 

With  these  initial  values  the  system  (25)  is  integrated  using  M-point  fourth-order  Runge- 
Kutta  integration.  The  general  solution  is  a  linear  combination  of  these  two  solutions; 

^  =  +  (30) 

where  are  chosen  to  satisfy  the  two  remaining  boundary  conditions  at  c  =  0.  Since 
the  boundary  conditions  at  2:  =  0  are  homogeneous,  this  requires  that  the  determinant  of  a 
certain  2x2  matrix,  whose  elements  are  the  appropriate  elements  of  the  vectors  —  1,2). 

vanish. 

However,  since  the  eigenvalue  problem  is  stiff,  rounding  errors  cause  the  base  solutions 
to  lose  their  linear  independence  fairly  quickly  during  the  integration.  To  circumvent 
this  difficulty,  the  base  vectors  are  orthonormalized  every  ten  integration  steps,  say  (a  total 
of  (2A//10)  —  1  orthonormalizations  since  no  orthonormalization  is  performed  at  the  final 
integration  step).  Thus,  every  ten  integration  steps,  is  replaced  by  where 

?'"’ =  II  (31) 

with  |[  ]|  in  the  denominator  denoting  the  modulus  of  the  bracketed  numerator.  In  satisfying 
the  boundary  conditions  at  7/1  =  0,  the  determinant  of  the  corresponding  2  x  2  matrix,  whose 
elements  are  now  the  appropriate  elements  of  ^  ,  (j  =  1,2)  is  required  to  vanish.  If  the 
eigenfunctions  are  needed,  it  is  necessary  to  keep  an  account  of  all  the  orthnormalizations. 
However,  to  determine  the  eigenvalues  it  is  only  necessary  to  solve  iteratively  the  non-linear 
equation, 

f{a,^,R,c)  =  0  (32) 

corresponding  to  the  vanishing  determinant  at  t/j  =  0.  This  is  done  here  using  Muller's 
method. 

4.  Results 

We  have  restricted  consideration  to  solutions  corresponding  to  the  most  unstable  dis¬ 
turbances.  The  neutral  curves,  which  represent  the  projections  onto  the  a  —  R  plane  of 
the  contours  of  the  neutral  surface,  given  by  ^  =  constant,  corresponding  to  three  values 
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of  the  spanwise  stations  Zq  {z*  =  lR~^zo  with  x*  =  1)  are  sliown  in  figure  2(a)-  (c).  Con¬ 
tours  are  plotted  for  values  of  /i  in  the  range  0  <  <  -Cr-  where  :i,.r  is  that  value  of  ,.i 
which  corresponds  to  the  critical  Reynolds  nuinl)er,  As  is  expected,  proximity  to  the 

corner  has  a  destabilizing  effect.  In  the  cases  shown,  instability  to  thre('-dimensional  distur¬ 
bances  in  the  above  range  of  spanwise  wavenumber  occurs  at  a  lower  Reynolds  number  than 
it  does  for  a  two-dimensional  disturbance,  this  feature  being  more  pronounced  at  stations 
closer  to  the  corner  than  those  away  from  it.  Proximity  to  the  corner  also  has  a  destabi¬ 
lizing  effect  on  the  two-dimensional  mode  of  disturbance  {R  =  0).  It  may  be  noteil  that  at 
.^0  =  6,  terms  neglected  in  the  basic  flow  in  (14)  may  become  important  and  inclusion  of 
these  terms  will  modify  the  stability  characteristics  at  this  location.  However,  we  expert 
that  the  characteristics  depicted  in  figure  2(c)  will  still  be  approximately  correct.  For  values 
of  2o  <  6  approximately,  we  expect  that  it  will  be  necessary  to  consider  the  exact,  fully 
three-dimensional  basic,  flow. 

It  is  interesting  that  the  instability  is  driven  by  the  outer  boundary  conditions,  the  most 
unstable  disturbance  having  a  finite  value  of  R.  This  is  consistent  with  Dhanak  (1981)  where 
it  was  shown  that  spanwise  waviness  in  one  of  the  walls  of  a  channel  has  a  destabilizing  effect 
on  the  plane  Poiseuille  flow  through  the  channel.  The  effect  is  much  more  pronounced  in 
the  present  case  in  view  of  the  presence  of  the  other  wall. 

The  neutral  curves  for  R  =  Rcr  associated  with  four  spanwise  stations  are  shown  in  figure 
3  for  comparison.  The  neutral  curve  for  a  Blasius  boundary  layer  [Rev  =  0)  is  also  included 
in  the  figure.  (Note,  however,  that  in  this  case,  contrary  to  normal  practice,  the  0(R~^) 
normal  component  of  the  basic  flow  is  retained  in  the  associated  Orr-Sommerfeld  equation, 
consistent  with  (20);  this  has  the  effect  of  shifting  the  neutral  curve  for  this  case  slightly 
towards  higher  Reynolds  number.)  The  critical  Reynolds  number  decreases  and  the  range  of 
the  unstable  streamwise  and  spanwise  wavenumbers  increases  as  the  corner  is  approached. 
The  values  of  the  critical  parameters  are  given  in  Table  1.  The  values  of  the  critical  Reynolds 
number  of  60  for  Zq  =  6  and  1  fO  for  ~o  =  7.5  are  consistent  with  the  experimental  prediction 
of  transitional  Reynolds  number  of  /?  =  100  by  Zamir  (1981).  The  value  of  R^t  increases 
while  acr  decreases,  so  that  the  critical  wave-angle  (cr  =  tan“' (/irr/arr)  increases,  as  the 
corner  is  approached; ter  =  48®  for  2o  =  6  compared  with  =  5°  for  =  ‘20.  As  can  be 
seen  from  figure  3,  the  neutral  curves  for  finite  values  of  Cq  approach  the  neutral  curve  for 
the  Blasius  boundary  layer  as  2:0  becomes  large. 

It  is  expected  that  the  most  unstable  mode  of  disturbance  will  dominate  the  initial 
development  of  the  instability.  The  growth  rate  of  the  most  unstable  disturbance  is  plotted 
as  a  function  of  the  Reynolds  number  for  each  spanwise  station  Zq  in  figure  4.  At  ro  =  6  the 
growth  rate  increa.ses  from  o-c,  =  0  at  /?  =  54  to  nc,  =  1 1.3  x  10“  ^  at  R  =  600.  At  R  =  600 
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the  growth  rate  of  the  disturbance  is  over  six  times  greater  at  ry  =  6  than  it  is  at  Cy  =  -0. 

The  variation  of  the  wave-angle  t  =  ta?i~'(/://A)  corresponding  to  the  most  unstable 
mode  of  disturbance  with  Reynolds  number  is  shown  in  figure  5.  Tlie  value  of  this  wave- 
angle  increases  as  the  corner  is  approached;  at  H  =  bOO,  (  =  11"  for  Cy  =  b  compared  with 
f  =  5'’  at  Zq  =  20.  f  decreases  at  first  as  H  is  rlecreasetl  from  its  value  of  bOO,  l)ut  then 
increases  somewhat  as  the  critical  Reynolds  number  is  ai)proached,  the  final  increase  being 
most  significant  for  the  case  ry  =  b;  c  =  48"  at  the  critical  Reynolds  number  for  Zq  =  b. 

5.  Conclusions 

It  is  shown  that  the  instability  of  the  flow  in  the  vicinity  of  a  streamwise  right  angled 
corner  is  dominantly  driven  by  oblique  disturbances  which  are  anti-symmetric  in  the  bisector 
plane.  The  magnitude  of  the  cross  flow  is  found  to  be  too  small  to  be  a  significant  factor 
in  the  observed  early  transition  of  this  flow.  It  is  found  that  close  to  the  corner,  three- 
dimensional  disturbances  in  a  particular  range  of  spanwise  wavenumbers  are  unstable  at  a 
lower  Reynolds  number  than  two-dimensional  disturbances.  Close  to  the  corner,  the  critical 
Reynolds  number  has  a  value  of  around  R  ~  50  compared  with  R  ~  300  for  the  Blasius 
boundary  layer  far  away  from  the  corner.  This  is  consistent  with  the  experimental  prediction 
(Zamir,  1981)  of  transitional  Reynolds  number  of  /?  =  100  for  this  flow.  At  /?  =  600 
the  growth  rate  of  the  most  unstable  disturbance  is  over  six  times  greater  at  a  distance 

=  6R~^x’‘  than  it  is  at  =  20R~^x*  while  the  corresponding  wave-angle  is  over  nine 
times  greater. 
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Table  1;  Critical  values  of  the  stability  parameters.  The  values  in  the  last  row  are  for  a 
Blasius  boundary  layer  and  are  evaluated  using  the  corresponding  Orr-Sommerfeld  equation 
in  which  terms  propotional  to  F,  v'  are  retained. 
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Figure  2.  Neutral  curves,  represented  by  contours  of  the  span  wise  wavenumber^  =  constant 
in  the  a-/?  plane,  at  spanwise  locations  (a)  Zo  =  20,  (b)  zo  =  10,(c)  Zo  =  6  .  Contours 
are  shown  for /3  in  the  range  0  S  ^  SI  where  is  the  value  of^  corresponding  to  the 
critical  Reynolds  number.  Here  spanwise  distance  z'  =  x’R~'zo, . 
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Figure  2(b). 
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to  a  Blasius  boundary  layer  is  included  for  companson. 
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Figure  4.  Growth-rate  of  most  unstable  disturbance  for  various  spanwise  locations  as  a 
function  of  Reynolds  number/?.  The  fixed  value  of/3  corresponding  to  maximum  growth 
rate  in  each  case  is  shown  on  the  associated  curve. 
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Figure  5.  Wave-angle  for  the  most  unstable  disturbiurce  as  a  function  of  the  Reynolds 
number  at  various  spanwi.se  locationsTo  .  The  corresponding  wave-angle  is  zero  for  a 
B  Iasi  us  boundary  layer  flow. 
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